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ABSTRACT 

We present high spatial resolution 21 cm HI observations of EA01A and EA01B, a pair of interacting 
post-starburst, or E+A, galaxies at z = 0.0746. Based on optical HST/WFPC2 images, both galaxies 
are known to display dist urbed morphologie s. They also appear to be linked by a bridge of stars. 
Previous HI observations (jChang et a l. 2001) had already uncovered sizable quantities of neutral gas 
in or near these galaxies but they lacked the spatial resolution to locate the gas with any precision 
within this galactic binary system. We have analysed deep, high resolution archival VLA observations 
of the couple. We find evidence for three gaseous tidal tails; one connected to EA01A and two 
emanating from EA01B. These findings confirm, independently from the optical imaging, that (i) 
EA01A and EA01B are actively interacting, and that, as a consequence, the starbursts that occurred 
in these galaxies were most likely triggered by this interaction, and that (ii) 6.6 ± 0.9 x 10 9 M Q of 
neutral gas are still present in the immediate vicinity of the optical bodies of both galaxies. The H I 
column density is lowest at the optical positions of the galaxies, suggesting that most of the neutral 
gas that is visible in our maps is associated with the tidal arms and not with the galaxies themselves. 
This might provide an explanation for the apparent lack of ongoing star formation in these galaxies. 
Subject headings: galaxies: evolution, galaxies: elliptical, galaxies: ISM, galaxies: fundamental pa- 
rameters 



1. INTRODUCTION 

Post-starburst galaxies (or PSGs, or E+A galaxies, 
or k+a/a+k galaxies) are characterized by their opti- 
cal spectra, independent of any morphological or pho- 
tometric considerations. PSGs contain a sizable pop- 
ulation of very young stars, producing very prominent 
Balmer lines, but lack ongoing star formation, hence 
the absence of detectable emission lines. This suggests 
that PSGs are indeed observed very shortly, within less 
than 1 Gyr, afte r the end of a vigorous, abruptly trun- 
cated starburst (IDressler fc Gunnl 119831 : iDressler et al.1 
119991 iPoggianti et al.lfl999h . Today, they constitute only 
a s mall fraction of t he ga laxy cluster population (< 
1%. iFabricant et al. I ([1991)); at intermediate-redshifts, 
howeve r, they constituted a substantial cluster popu- 
lation (|Belloni et a l. 1995). Although the trigger and 
the abrupt end of the starburst is still not fully un- 
derstood, photometry of PSGs shows evidence for dis- 
turbed morphologies, e.g. tidal tails, suggesting that 
in many cases the trigger is most like ly associated to 
a gal a xy-galaxy merger or interaction | Zabludof f et al 
1996; Yang et al. 2004; Blake et al.l 120041 : iTran et al 
20031: iGotol I2005t lYang et all I2008T ). Using numerical 
simulations, iBekki et al.l (|2005[ ) has shown that PSGs 
can be formed via a major merger of two gas-rich spi- 
ral galaxies. In these merger simulations, a starburst is 
triggered that consumes the available gas within a times- 
pan of roughly 1 Gyr, ending the starburst abruptly. Re- 
cent H I observations have uncovered large quantities of 
neutral hydrogen gas in a signif icant fraction of PSGs 
(IChang et al. 1120011 : Huvle et al.ll200J iHelmboldt 1120071: 
iBuvle et al.ll2008l ). This opens up the possibility of inves- 
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tigating the hypothesized merger origin of many PSGs. 

Indeed, high-resolution interferometric radio observa- 
tions of the neutral gas would provide valuable additional 
information. On the one hand, H I morphologies and ve- 
locity fields, if they are distorted, could corroborate the 
merger hypothesis, and, moreover, could reveal whether 
at least part of the neutral gas is still gravitationally 
bound and will remain available for future star forma- 
tion. We selected the PSG binary EA01A/B as an inter- 
esting target for our study. Both galaxies have a PSG- 
type optical spectrum and have the same recession ve- 
locity. EA01A i s the b luest galaxy of the PSG sample of 
IZabludoff et al.1 (|1996l ). suggesting it is also the youngest 
sample member. HST/WFPC2 im ages in the F435W 
and F702W bands (|Yang et al. 2004) provide strong evi- 
dence that the couple is interacting. EA01A, positioned 
about 11" east of EA01B, contains many stellar clus- 
ters that are very conspicuous compared with the overall 
rather low surface brightness of this galaxy. This diffuse 
appearance suggests it is on the verge of disintegrating 
due to the injected orbital energy. EA01B, on the other 
hand, is a bulge dominated early-type spiral galaxy; it 
sports strongly asymmetric stellar arms. Both galaxies 
appear to be connected by a stel lar bridge. Previous HI 
observations (|Chang et al. l|200i] ) had already uncovered 
sizable quantities of neutral gas in or near these galaxies 
but they lacked the spatial resolution to locate the gas 
with any precision within this galactic binary system. 

In section [3J we describe our data reduction and 
analysis of deep archival VLA H I observations of the 
EA01A/B galaxy binary. The results are presented and 
discussed in section [3] and summarized in section 2] 

2. OBSERVATIONS AND DATA REDUCTION 

We retrieved deep, high-resolution H I 21 cm obser- 
vations of EA01A/B from the science archive of the 
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Fig. 1. — Channel maps overplotted on an optical HST/WFPC2 F702W image of EA01A/B. The contour levels are 2a, 3a,. . .with 
=0.2 mjy beam -1 . The synthesized beam is depicted in the bottom-right corner of the first map (top left). 



TABLE 1 
Data characteristics 



channel width (km s -1 ) 

synthesized beam (FWHM) 

la noise level per channel 

flux density (mjy beam -1 ) 
brightness temperature (K beam -1 ) 
column density (xlO 19 cm -2 beam -1 ) 

integrated la noise level 

flux (Jy beam - km s -1 ) 

brightness temperature (K beam -1 km s — x ) 

column density (xlO 20 cm -2 beam -1 ) 



11.1 

8.87"x6.32" 

0.2 
2.1 
4.3 

0.01 
71.1 
1.9 



National Radio Astronomy Observatory (NRAO), made 
with the Very Large Array (VLA) in New Mexico (USA) 
(project number AM0678). The observations were made 
with the B configuration during the nights of March 
11, 12, 14, 16 and 18, 2001. At a central frequency 
of 1321.76 MHz (L-band) and with a bandwidth of 



3.125 MHz divided over 63 channels, the H I in EA01A/B 
was mapped with a velocity resolution of 11.1 kms -1 . 
Alternating observations of the phase calibrator 1127- 
145 with a 4 minute exposure time and EA01A/B with 
an exposure time of 25 minutes were made, followed at 
the end of each night by a 2.5 minute exposure of the flux 
calibrator 3C286. This resulted in a total on-source inte- 
gration time of 32.6 h. Standard flagging and calibration 
of the u — v data was performed with the Astronomical 
Image Processing Software (AIPS). The continuum was 
subtracted by making a linear fit to the visibilities over 
the line-free channels that were not affected by the edge 
effects of the band. We created the final datacube using 
natural weighting to optimize the sensitivity. This yields 
a beam-size of 8.87" x 6.32" or 12.3 hr x kpcx8.7 fr" 1 kpc 
at the angular-size distance of EAl (D^ = 285 h^ 1 Mpc, 
with Ho = 75 h kms" 1 Mpc -1 ). The rms noise per chan- 
nel is 0.2 mjy beam -1 , corresponding to a H I column 
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Fig. 2. — The observed H I flux integrated map (left panel) from Chang et al. (2001) (white contours) and from this paper (black 
contours) and the H I flux integrated map of a model that roughly reproduces the most conspicuous features of the channel maps (right 
panel). The synthesized beam is shown at the bottom left. The contour levels of the observed flux map correspond to 2cr, 3a,. . .with 
a = 0.01 Jy beam -1 kms -1 . The model channel maps are displayed in the bottom panels. 



density of 4.3 x 10 cm averaged over the beamwidth. 
A summary of the data characteristics can be found in 
Table [U 

3. DISCUSSION 

We examined the final datacube and detected 21 cm 
emission in 12 adjacent channels, shown in Fig. [1] The 
H I flux integrated map of EA01A/B is presented in the 
left panel of Fig. O All maps are overplotted onto a 
HST/WFPC2 F702W image of the EA01A/B pair. It 
is clear from these maps that the H I emission spatially 
coincides with the optical position of EA01A/B and is 
centered around the heliocentric velocity of the previous 
H I detection and the optical recession velocity of the 
pair (22327.9 - 22450.0 kms" 1 ). 

The H I profile of EA01A/B, obtained by summing 
the flux within a box covering the radio emission of the 
galaxy is plotted in Fig. [3l We fitted a Gaussian to 
the H I spectrum in order to derive the total 21 cm 
line flux. We generated a large library of mock H I 
spectra using this Gaussian model and added noise to 
produce synthetic data of the same quality as the orig- 
inal spectrum. We then fitted Gaussian profiles to the 
mock data in order to estimate the la uncertainties on 
all relevant parameters. We find a total integrated flux 
of S u = 0.26 ± 0.03 Jy kms -1 . This is equivalent to 
an H I mass of Mhi = 6.6 ± 0.9 x 10 9 M Q , adopting a 
luminosity distance Di=330 Mpc. The central velocity 
of the Gaussian is 22357 ±9 kms" 1 and its FWHM is 
133 ± 20 kms -1 . We also re-redu ced and analy s ed the 
low-resolution H I observations of IChang et al. I (2001) 
(Fig. [2]). Based on these data, we find a total integrated 
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Fig. 3. — The H I profile of EA01A/B, obtained by summing 
the flux within a box around the flux integrated radio emission of 
the galaxy. Rms noise levels are indicated by horizontal lines. A 
Gaussian was fitted to the spectra to measure the velocity widths. 



flux of S v = 0.29 ± 0.08 Jy kms 1 , corresponding to an 
H I mass of 7.5 ± 2.1 x 10 9 M©. The central velocity is 
22371 ± 30 kms -1 and the FWHM is 191 ± 88 kms -1 . 
This is in agreement wit h the high-res o lution VLA data 
presented here and with lChang et alTl (|2001[ ). 
A plausible interpretation of the H I flux integrated 
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map and the channel maps is that we are seeing one re- 
ceding gaseous tidal arm emanating from EA01A and two 
tidal arms connected to EA01B, with the east arm ap- 
proaching us and the west arm moving away from us. Lit- 
tle of the detected gas appears to be directly connected 
with the galaxies themselves, which show up as minima 
in the H I column density. This absence of dense gas may 
be related to the fact that these galaxies are in a post- 
starburst phase. The EA01A arm runs in a south-eastern 
direction and can be traced through 3 consecutive chan- 
nel maps from 22372.3 km s -1 to 22394.5 km s -1 . The 
west EA01B "arm" shows up as emission in 3 consecutive 
channel maps from 22372.3 km s" 1 to 22394.5 km s _1 . 
Following its position through adjacent channel maps 
with recession velocities below rs 22370 km s -1 , the east 
arm connected with EA01B appears to first shift in the 
direction of EA01 A and then to return a position slightly 
north of EA01B before disappearing in the noise. At that 
point, the gas is moving towards us with a radial velocity 
of about 40 km s _1 , relative to the systemic velocity of 
EA01A/B. The curvature of the northern arm and the 
fact that it is approaching us agrees with the curvature 
of the stellar arms. 

Since the data quality clearly does not allow for a de- 
tailed or quantitative modeling of the gas distribution 
and kinematics, we constructed a very simplified model 
for the H I observations based on the interpretation out- 
lined above. Two gaseous arms were constructed to em- 
anate from EA01B, following the position and curvature 
of the stellar arms. If the ellipticity e = b/a rs 0.5 of the 
outer isophotes of EA01B is indicative of its inclination, 
then this early-type disk is inclined by about 30° from an 
edge-on view. We assumed both arms to lie in a disk with 
the same inclination. A third gas arm was connected 
with EA01A, curving in a south-eastern direction. The 
H I flux density of these model tidal arms was assumed 
to decline exponentially with radius and they were given 
solid-body rotation velocity fields in such a way that the 
mock channel maps and the flux integrated map derived 
from this model provided an acceptable match to the ob- 
served maps, taking into account beam convolution and 
added noise. Despite its simplicity, this model can give 



a fair account of the observations. The model flux map 
and channel maps are presented in Fig. [2] and can be 
compared with Fig. [TJ 

With a 2cr density limit of 4x 10 20 cm~ 2 beam -1 , aver- 
aged over a 12.3 hT l kpcx8.7 h~ l kpc beam, we cannot 
expect to trace the tidal arms very far out. As a com- 
parison, under these conditions it would be impossible 
to trace the southern arm emanating from the interact- 
ing pair NGC40 38 /4039 ( "The Anten nae" ) further out 
than ~ 20 kpc (|Hibbard et al. 1120011 ). At the distance 
of EA01A/B, this corresponds to ~ 15", less than two 
beam diameters. At that distance, the velocity of the 
NGC4038/4039 tidal arm is less than - 50 km s"\ in 
rough agreement with the largest velocities measured in 
the EA01A/B system. 

4. CONCLUSIONS AND SUMMARY 

EA01A/B is a close pair of post-starburst (E+A) 
galaxies, surrounded by some 7 x 10 9 M of neutral gas. 
Most of this gas resides in what appear to be three tidal 
arms, two of which are connected with EA01B. Together 
with optical HST/WFPC2 images, these observations 
show that EA01A and EA01B are actively interacting. 
The galaxies themselves show up as minima in the H I 
column density. This lack of galaxy-bound dense neu- 
tral gas is most likely connected with the fact that these 
galaxies are in a post-starburst phase. 
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